Phonon structure in point-contact spectra of MgB 2 
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In strong-coupling superconductors with a short electron mean free path the self-energy effects in 
the superconducting order parameter play a major role in the phonon manifestation of the point- 
contact spectra at above-gap energies. We compare the expressions for the nonlinear conductivity 
of tunnel, ballistic, and diffusive point-contacts and show that these expression are similar and 
correspond to the measurements of the phonon structure in the point-contact spectra for the ir- 
band ofMgB 2 . 
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I. INTRODUCTION 

The inelastic point-contact spectra of phonons in met- 
als are based on expansion of the nonlinear I — V char- 
acteristic in terms proportional to d/l e , where d is the 
characteristic size of metallic nanoconstriction connect- 
ing two bulk metal half-spaces and l e is the inelastic elec- 
tron mean free pathA^. 

If d is smaller than the electron mean free path I = 
l^ 1 + l^ 1 , k is the elastic mean free path), the electrons 
flow through the contact ballistically. In opposite cases 
(d 3> I) they perform a diffusive motion, but if the energy 
relaxation length A £ m ■s/TJe remains larger than d, their 
movement inside the constriction is conservative, and the 
inelastic point-contact spectroscopy is still possible^. 

In the superconducting state a new characteristic 
length appears: the size of the Cooper pair or the co- 
herence length £, and the additional nonlinear fea- 
tures emerge being proportional to the expansion of 
the so called "excess" current in terms proportional to 
d/£. In all the cases considered, the size of the contact 
is smaller than the corresponding characteristic length 
scale: £, l e , A e etc. 

The zero-approximation in expansion of /(V) in d/l e 
is the ohmic current considered by Sharvin^. From the 
Sharvin's formula the characteristic size d of the con- 
tact can be determined in the ballistic limit. The sec- 
ond derivative of the first approximation in expansion 
of I(V) in djl e is directly proportional to the spec- 
tral function of electron-phonon interaction (PC EPI) 
gpc(uj) = ap c (uj) F (uj) of the specific point-contact 
transport both in the normal and in the superconducting 
statesiaSi This term is the basis of the canonical inelas- 
tic point-contact spectroscopy (PCS). Here, ap c (uj) is 
the average electron-phonon matrix element taking into 
account the kinematic restriction imposed by contact ge- 
ometry and F (uj) is the phonon density of states. 

In diffusive point contacts (d^S> U) the role of the scale, 
where the backscattering inelastic processes become es- 
sential, turns from d in and in the case when li/l e <C 
1 the corresponding contribution becomes negligible^. 
This property will be essential when we consider the 
phonon structure in point-contact spectra of dirty MgB2 



contacts in the c-direction. 

In the superconducting state, due to the Andreev re- 
flections of quasiparticles at the N — S boundary, the 
excess current appears both in S — c — N and S — c — S 
contacts (c stands for constriction). Its magnitude is the 
constant proportional to A for the weak-coupling BCS 
approximation at eV 3> A (A is the superconducting en- 
ergy gap)^. In the strong-coupling model of EPI the 
energy gap depends on the energy of quasiparticles due 
to phonon retardation effects reflecting the characteris- 
tic features of the phonon spectrum at singular van-Hove 
energies. Hence, even in the zeroth order the nonlinear- 
ities appear in 7(V)iS. For the issue of this study it is 
most important that these nonlinearities are independent 
of the elastic mean free path (i.e., purity of the contact). 
In this connection one should keep in mind that for the 
simple method of a gentle touch between two sharp edges 
of the electrodes it is not possible to control the elastic 
mean free path in the constriction. In what follows we 
show that in spite of this restriction, the satisfactorily 
reproducible results are obtained not only for the shapes 
of phonon structure but also for their amplitude. 

To conclude this section, we mention that the use of 
the first approximation in expansion of excess current in 
powers of d/£ allows the nonequilibrium effects in the 
vicinity of the constriction to be developed, which are 
out of the scope in this paper—. 



II. THEORETICAL BACKGROUND OF 
ZEROTH-ORDER ELASTIC COMPONENT IN 
EXCESS CURRENT 



In this section we review what has recently become 
known about the phonon structure in the point-contact 
spectra—. In order to do this, we should refer to the 
works where the excess current in point contact is con- 
sidered in terms of the strong-coupling theory. That was 
first done in the paper of Omcl'yanchuk, Beloborod'ko 
and Kulik for ballistic S — c — N point contact. The 
first derivative of /(y)-characteristic at T = has the 
following form in this case: 
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where Rn is the resistance of the contact in the normal 
state. This expression implies d <C £ to be fulfilled (£ is 
the coherence length in pure metal). 

For diffusive point contact with the restriction: li <C 
d <C v^iCi the same quantity takes the formM: 
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We would like to compare these dependences with the 
expression for the conductance of the tunnel junction due 
to the self-energy effect in the superconducting energy 
gap: 



Rn[ 1v 



= 3? 



yje 1 - A 2 (e)' 



(3) 



In the above mentioned formulae A (e) is the complex 
superconducting order parameter 



A(e)=A 1 (e)+iA 2 (e) , 



(4) 



where Ai (e) and A 2 (s) are the 3? (e) and 3 A (e), respec- 
tively. It would be very instructive to express Eqs. 0J, 
(0, J21 as functions of real and imaginary components. 
Thus, using the auxiliary expression 



A (e) = •^(e 2 -A 2 + A 2 ) 2 +4A 1 A 2 
the results are as follows: 
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where "plus" and "minus" in the last line stand for 
AiA 2 < and AiA 2 > 0, respectively. 

Finally, for tunnel junction the same quantity reads: 
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Expressions 0, ©, © [and ©, Q, ©] seem quite 
different but in reality they are very similar. This be- 
comes evident if we compare their asymptotic expressions 
for the limit £> A. For three extreme regimes of current 
flow in aforementioned point contacts the corresponding 
equations take the form: 

for a ballistic point contact: 
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for a diffusive point contact: 
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and for a tunnel junction: 
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In all these equations e = eV and T = are implied. 

Since A 2 (e) becomes comparable to Ai (e) only at the 
highest phonon energy, almost in the whole phonon en- 
ergy range these three extreme regimes of current flow 
differ only little, mostly in amplitudes, which decrease 
in the tunnel— >diffusive—>ballistic contact row like their 
prefactors \ — > | — > |. This derivation is of a key 
importance to understanding the quite similar phonon 
point-contact spectra for the different contacts randomly 
accessible in the experiment. 

Let us illustrate the above theoretical formulae by 
an example of superconducting Pb frequently used for 
these purposes in tunneling spectroscopy of the electron- 
phonon interaction. This metal has the relatively strong 
EPI and a simple phonon spectral function, consisting 
of two acoustic peaks: transverse (TA) and longitudi- 
nal (LA) (Fig. 1). We take the tabulated functions: 
a 2 (e)F(e), 3?A (e) and 3A(e) from Rowell-McMillan- 
Dynes/ preprint 13 and plot both the EPI spectral function 
and the second derivative characteristic — {RnM?) 
proportional to the experimentally recorded second har- 
monic of modulation voltage, V2(eV) (see the text be- 
low). One can immediately see that three point-contact 
spectra (tunnel, diffusive and ballistic) are very similar 

"5RA(e) 



differing only in the constant prefactor at — (see 

above). The differences in shape occur only at the sec- 
ond peak of the EPI spectrum, where Q : A(e) becomes 
essential. For the transverse peak the positions of max- 
ima of point-contact spectra V 2 (eV) coincide with the 
TA-maximum of EPI function shifted by A to the right 
( we remind that throughout this paper S — N contact is 
considered). On the other hand, the maxima of V 2 (el/) 
functions at LA peak are slightly displaced with respect 
to the EPI function and their shapes change, which be- 
comes more pronouced the further we deviate from the 
tunnel junction in the row mentioned above. 
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FIG. 1: The second derivative of I — V characteristics for 
superconducting lead S — c — N point contact calculated 
by Eqs. ©, Qi ® f° r the tunnel (solid curve), diffusive 
(dash-dotted curve) and ballistic (dotted curve) regimes of 
current flow. 5RA (e) and SsA (e) are taken from the Rowell- 
McMillan-Dynes 1 - 3 tabulation. Squared symbols stand for 
Rowell-McMillan-Dynes— EPI spectral function shifted by 
A = 1.4 meV to the right on the voltage scale. Its ordinate 
scale is shown on the right side. 



Comparing the exact spectra (Eqs. JBJ, 0, ©) with 
the asymptotic ones (Eqs. ©, 11UI) . (|llfl ). we see that for 
Pb the differences are only minor (the example for diffu- 
sive point-contact is shown in Fig. 2). It goes without 
saying that the differences for materials with greater ra- 
tios of phonon energy to the energy gap, hu> p h/ A, will 
be very small. For example, in Sn or MgB2 , where 
hujph/A. 3> 1 the differences between exact and asymp- 
totic curves in the phonon energy range become indistin- 
guishable. 

Experimentally, the expected shift of the Pb phonon 
spectra to the right was first noticed in Ref. Unfor- 
tunately, in this paper (and also in the review paper—) 
the experimental curves were compared with theoreti- 
cal ones without allowing for the shift of the latter by 
A to the right. In Fig. 3 this drawback is corrected 
and the important consequences are drawn for the bal- 
listic point contact. The upper panel of Fig. 3 dis- 
plays two experimentally recorded dependences: i) the 
second derivative of the I_V characteristic in the nor- 
mal state (superconductivity is destroyed by the mag- 
netic field) and ii) the same dependence in the super- 
conducting state. The BCS dependence in the supercon- 
ducting state is shown for comparison. The theoretically 
expected point-contact spectrum for ballistic contact is 
shown in the bottom panel of Fig. 3. It is compared 
with the experimental PC EPI spectrum in the normal 
state shifted by Apt, = 1.4 meV . As expected, the po- 
sition of the calculated TA-maximum coincides with the 
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FIG. 2: Comparison of the exact Q and asymptotic HUH 
calculated curves for lead in the S — c — N point contact. 

TA-maximum in the normal-state spectrum and the cor- 
responding feature (zero of y-coordinate) coincides with 
the normal-state LA-maximum (see Fig. 1). 

In contrast, the phonon features on the experimental 
curve in the superconducting state shown in the upper 
panel are noticeably displaced to the left from the ex- 
pected calculated positions (shown in the lower panel) 
towards the non-shifted phonon maxima. This displace- 
ment is the result of mixing at least of two different mech- 
anisms of revealing the phonon features: namely, the in- 
elastic backscattering current, which has the same posi- 
tions both in the superconducting and the normal states 
(the latter is shown as an experimental PC EPI curve in 
upper panel of Fig. 3), and the theoretical zeroth-order 
approximation of excess current shown as a theoretical 
curve in the lower panel. Such mixing seriously impedes 
interpretation of PC spectra of a ballistic strong-coupling 
superconductor. Fortunately, as mentioned above, in dif- 
fusive point contact the inelastic backscattering current 
is negligible and the pure self-energy effect on excess cur- 
rent becomes clearly visible. The latter situation holds 
in MgB2 films measured along the c-direction as well. 



III. SELF-ENERGY EFFECT IN c-AXIS 
ORIENTED MGB 2 FILMS 

For two-band superconductor MgB2 (see the reviewi^) 
the SKA (uj) and OA (uj) for tt and u-band are depicted in 
the upper and lower panel of Fig. 4i£, respectively. We 
are interested mostly in the 7r-band dependences, since 
for thin films oriented along the c-axis only the 7r— band 
part of the Fermi surface has nonzero Fermi velocity for 
charge carriers (holes) along the contact axis, playing the 
main role in the contact current. For the cr-band, which 
is two-dimensional, only holes with the velocity in the 
ab-plane exist. The driving force for superconductivity is 
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FIG. 3: The Pb-Ru experimental point-contact spectra in the 
normal (squared symbols) and superconducting (triangles-up) 
states (upper panel)—. The calculated spectrum for a ballistic 
point contact is shown in the lower panel (triangles-down) 
together with the experimental PC EPI spectrum shifted by 
A = 1.4 meV to the right. The theoretical BCS dependence 
is shown for comparison (thin dotted curve). 



the very strong EPI of these holes with the i^g-phonon 
mode with the energy of about 70 meV. This mode is seen 
as the strong peak in SftAo- in the lower panel of Fig. 4 15 . 
Note, that the ordinate scale for the latter panel is about 
three times higher than that for the 7r-band. For the 
7T— band all the phonons in the energy range ~ 30 ± 100 
meV take part in the SRA^ dependence, but with a lower 
magnitude as compared with the er-band. Up to about 
80 meV SA (o>) is small, hence, only 5RA (ui) is important 
in Eqs. ©, EJ, flTTJ). Thus, we expect that up to ~ 80 
meV the point-contact spectra differ onlj in amplitude 
just like the prefactors in the [9?A (e) /e] term. This is 
evident from Fig. 5. 

We can see that up to the energies ~ 80 -j- 90 meV 
the shapes of the calculated spectra are similar (Fig. 
5). The discrepancy is observed only at higher energies, 
where 9A (cu) becomes essential. In the range of biases 
80-^120 meV the variations of point-contact spectra are 
quite appreciable and mirror in the experimental curves. 
Although the A (w) value is of the order of A max at bi- 
ases higher than 120 meV, the phonon features on the 
PC spectra are small enough, due to the factor 1/e 2 (u>). 



FIG. 4: Real and imaginary parts of A (e)-dependences for 
7r- and cr-bands (upper and lower panels, respectively) of 
MgB;A 



When we compare theoretical and experimental curves, 
we consider only the biases up to 120 meV, since for 
higher biases the experimental results are irreproduciblc, 
due to nonequilibrium suppression of the superconduct- 
ing state in the contact. 

Comparison of calculated and experimental curves is 
shown in Fig. 6. For the ordinate scale we use the second 
harmonic voltage V% measured in experiment, which is 
related to the second derivative of the I—V characteristic 
as follows: 



V 2 (V) 



2V2dV 



(12) 



where V\ is the modulation voltage, which is convention- 
ally taken equal to 3 meV (close to the value used in the 
experiment). 

In the upper panel of Fig. 6 three different point- 
contact spectra are shown. They have similar shape up 
to 100 meV with a slight increase of average amplitudes 
in the row Rn = 49-^36-^80 f2. In the same sequence 
the so-called zero-bias maximum in the dV/dI{V) char- 
acteristic increases (not shown) , being Rq/Rn < 1, ~ 1 
and > 1, respectively (i?o is the zero-bias resistance in 
the superconducting state), for the above mentioned row. 
One may speculate that in the same direction the con- 
tact changes from ballistic towards the tunnel regime 22 . 
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FIG. 5: Calculated point-contact spectra for n- and er-bands 
(upper and lower panels, respectively) of MgBz— for ballistic 
(dashed curve), diffusive (thick solid curve) and tunnel (dot- 
ted curve) regimes of current flow. Real and imaginary parts 
of A (e) are taken from Fig. 4. 



The lower panel of Fig. 6 shows that in the same row the 
average amplitude of phonon structure increases in the 
calculated curves. 

The positions of the maxima of the experimental 
phonon features roughly correspond to the expected ones 
taking into account that the orientation of experimen- 
tal junction is approximately along the c-axis, since the 
spreading behavior of current in direct conductivity point 
contact embraces a wide solid angle near the contact axis. 
The alignment of our contacts along the c-axis is inferred 
indirectly by measuring the Andreev reflection spectra 
in the energy gap region and by observation of only the 
small gap A w 2.4 meV characteristic of the 7r-band. 

Special attention should be paid to the amplitude of 
the experimental curves, which roughly equals that pre- 
dicted theoretically. In inelastic spectroscopy the ampli- 
tude of the EPI spectrum is an order of magnitude lower 
than expected one (see Fig. 4 in Ref. and the discus- 
sion cited therein). This discrepancy may be explained 
either by the diffusive regime of current flow with and un- 
known mean free path ij, or by the specific PC-transport 
character of the EPI function obtained from the inelastic 
backscattering spectrum. 

Comparing theoretical and experimental spectra, one 
can infer that all the phonons are essential in the EPI 



FIG. 6: Comparison of experimental (upper panel) and calcu- 
lated (lower panel) phonon structures in point-contact spectra 
of MgB2 . The contact axis is oriented along the c-axis. The 
normal state resistances are given for each curve. In the ex- 
perimental panel the upper curve corresponds to the ordinate 
scale, and the other two are shifted down for clarity. T = 4.2 
K. In the theoretical (lower) panel three transport regimes are 
illustrated on the same scale as in the experimental graph. 
Here, the lower curve corresponds to the ordinate scale, the 
other two are shifted up for clarity. The 7r-band EPI function 
(dashed curve)— is shifted to higher voltages by A^ = 2.4 
meV. The modulation voltage is taken equal to 3 mV. 

function for the 7r-band. In the lower panel of Fig. 6 we 
plot the (a 2 F)^ function taken from Ref. El As ex- 
pected, the maxima in the EPI spectral function shifted 
by A^O) to the right coincide approximately with the 
maxima in the second harmonic dependences V2(V). 



IV. CONCLUSION 

The self-energy effect in the phonon feature of a su- 
perconducting point contact can be used, in principle, 
in the same way as the Rowell-McMillan program for 
determination of the EPI spectral function in tunneling 
spectroscopy of superconductors^. Two difficulties arise 
on this way. One is theoretical, since this program works 
well only for the one-band superconductor, and its ap- 
plication to the two-band case, like MgB2, encounters 
difficulties^. The other is experimental, since all other 
sources of J — V nonlinearities should be removed, and 
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especially, the nonequilibrium effects in superconductor 
should be excluded. 

The result obtained for superconductivity in MgB2 is 
the first direct proof of validity of the calculated EPI 
function- 7 , at least in the 7r-band. Unfortunately, this 
method encounters difficulty, while applying it to the a- 
band. In this case the very strong generation of nonequi- 
librium phonons cannot be efficiently excluded from the 
contact region and this destroys the superconductivity in 
the constriction. 

The question may arise whether the self-energy effects 
are important in the normal state. These are known to be 
smaller than the inelastic backscattering nonlinearities in 
the ballistic regime^. If we decrease the contact size d or 
the elastic mean free path U in order to make the inelas- 
tic contribution negligible, the latter parameters become 



comparable to the Fermi wave length of charge carriers 
and the strong nonlinearities connected with localization 
occur, which masks the desired phonon structure^. 
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current flow. This happens in many cases, where the same 
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